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We present a study of the pseudoscalar and vector meson form factors, calculated using the Fat- 
Link Irrelevant Clover (FLIC) action in the framework of Quenched Lattice QCD. Of particular 
interest is the determination of a negative quadrupole moment, indicating that the p meson is not 
spherically symmetric. 
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I. INTRODUCTION 

The important role that electromagnetic form factors 
play in our understanding of hadronic structure has been 
well documented for more than fifty years. The reason 
for their popularity is that they encode information about 
the shape of hadrons, and provide valuable insights into 
their internal structure in terms of quark and gluon de- 
grees of freedom. 

Most of the attention, both experimentally and theo- 
retically, has focused on the electromagnetic form factors 
of the nucleon (see Refs. 0, 0, S H H for recent reviews). 
The electromagnetic form factors of pseudoscalar mesons, 
especially the pion, being the lightest QCD bound state, 
have also been studied extensively [1, @, H, H EH EJ 
in lattice QCD. More recently, there is a renewed in- 
terest in calcula ting the pion form factor on the lattice 
El E Q El Ea 111 El- This is especially timely 
considering the new [19( and reanalysis of old [2fJ exper- 
imental data from JLab. 

The vector meson form factors, on the other hand, have 
received less attention (see Refs. [2l|, [22|, [23|, HH HH for 
recent work). Of particular interest is the quadrupole 
moment of the p meson, where theoretical determinations 
can disagree by as much as a factor of two [22j ■ We aim to 
resolve this issue by performing the first direct lattice cal- 
culation of the p- meson quadrupole form factor. Charge 
and magnetic form factors are also calculated and from 
these we extract the relevant static quantities, namely 
the mean square charge-radius and magnetic moment. 
We also analyse the dependence of light-quark contribu- 
tions to these form factors on their environment and con- 
trast these with a new calculation of the corresponding 
pseudoscalar-sector result. 

Our aim is to reveal the electromagnetic structure of 
vector mesons and to study to what extent the qualitative 
quark model picture is consistent with quenched lattice 
QCD. Interestingly, it has been shown in a lattice calcu- 
lation by Alexandrou et al. [26| that the distribution of 
charge in the vector meson is oblate, and therefore not 
consistent with the picture of a quark anti-quark in rela- 
tive S-wave. By calculating the vector meson quadrupole 



form factor we make a direct comparison with the find- 
ings of Ref. [H. 

For each observable we calculate the quark sector con- 
tributions separately. Using this additional information 
we examine the environmental sensitivity of the light- 
quark contributions to the pseudoscalar and vector me- 
son charge radii. We also evaluate the dominance of the 
light quark contributions to the K and K*. 

This paper builds on the preliminary work presented in 
Ref. [27|] . In Section Bl Al we introduce the theoretical for- 
malism of meson form factors, including the techniques 
required to extract them from a lattice calculation. Sec- 
tion 3 contains details of our lattice simulation, while 
in Section 4 we present and discuss our results for both 
pseudoscalar and vector mesons. Finally, in Section 5 we 
summarise our findings and discuss future work. 



II. THEORETICAL FORMALISM 

A. Meson form factors 

Meson form factors are extracted from matrix elements 
involving the vector (electromagnetic) current 



(M(p')\J a \M(p)} 



(1) 



where M(p) (M(p')) denotes a meson state with initial 
(final) momentum p (p'). The momentum transfer is 

For a pion, the matrix element in Eq. {T]) is described 
by a single form factor 



(■K{p')\J a \ir{p)) = 



[p a +p a ']F n (Q 2 ),(2) 



where Q 2 = —q 2 is the invariant momentum trans- 
fer, and the energy of the pion with momentum p is 
E-n{p) = \J m 2 + p 2 . The formula for the kaons are ex- 
actly analogous. The p- meson, on the other hand, is 
spin-1 and is described by three form factors [2§1, 

(p(p',s')\J a \p(p,s)) = 



2 



e'*(p', s') J Taa (p',p) e a {p, sX3) tors ' respectively, and 



2^E p (p)E p (p<) 
where e and e' are the initial and final polarisation vec 



J Ta °{p\p) = ~ Gi(Q ) g™ [p a +p al ] + G 2 {Q')[g a °q T - g aT q°\ - G 3 (Q 2 ) q°q 



a r P~+P" 



Gq(Q 2 ) = G^Q 2 ) - G 2 (Q 2 ) + (1 + V )G 3 (Q 2 ) (5) 



G M (Q 2 ) = G 2 (Q 2 ) 

G C {Q 2 ) = G 1 (Q 2 ) + ^r / G Q (Q 2 ). 



(6) 
(7) 



The charge q p , magnetic moment /_t p , and quadrupolc 
moment Q p are then extracted from Gc, Gm, and Gq, 
respectively at zero momentum transfer 



eGc(O) = q p 
eG A /(0) = 2m p p, p 
e Gq(0) = m 2 " 



1 p^p ■ 



(8) 
(9) 
(10) 



The formulae for the K* are exactly analogous. 



B. Meson form factors on the lattice 

The matrix elements in Eqs. ((2]) and ((3]) arc obtained 
from ratios of three-point and two-point correlation func- 



2m 2 



(4) 



r 



The covariant vertex functions 61,2,3 can be rewritten tions 
in terms of the Sachs charge, magnetic and quadrupolc 
form factors [28L \2M , 



R a (p',p) 



'(G«{p',pMM)) (G«(p,p'mm)) 



<G(p> M)><G{pM)> 
for pseudoscalar mesons, and 



(11) 



where m p is the mass of the vector-meson system calcu- R a v (p ,p) 
lated on the lattice and 7/ = Q 2 /Am 2 . 



I (Gg„(p' JMM)) (Gg M (p,p',f 2 ,ti)) 
<G flfl (t 2 ,p')> <G vv (t 2 ,p)> 



(12) 

for vector mesons. Note repeated indices are not summed 
over. 

The two-point correlation function for the pseudoscalar 
mesons is 

G{t 2 ,p) = ^e-^ 2 (f2|x(^ 2 )x t (0)|Q) . (13) 



Similarly for the vector mesons, 

G^(t 2 ,p) = ^e-^ (il\ Xp (x 2 ) X t(Q)\n) . (14) 



The three-point correlation function for the pseudoscalar meson is 

G a {t 2 ,h,p',$) = £ e -<?'-(* a -*0 e -*?-*i (0| x ( X2 )J°( Xl ) x l(0) |fi) . 

Similarly the three-point function for the vector meson is 

G%„(t2,ti,p',p) = Yl e-W'-^e-** 1 (fi| Xp(x 2 )J a (xi)xt(0) |0> • 



(15) 



(16) 



The Lorentz indices fi and v are only present for the 
vector mesons, while a is the index of the electromagnetic 
current. 

The ratios in Eqs. tfTTj) and (fT2"]l are constructed in such 
a way as to remove the time-dependence and constants 



r 



of normalisation from the correlation functions at large 
time separations, t\ and t 2 —t\. 

These ratios differ subtly from previous work (30| , in 
that we are explicitly enforcing the parity of the terms 
through the choice of momenta (p' ,p) and (p,p') vs (p 1 ,p) 



3 



and (— p, —p'). This requires two three-point propagators 
(with momentum-transfer q and —q) for each configura- 
tion. However with the well established technique of av- 
eraging over U and U* configurations [13, there is no 
additional cost. 



1. it -meson case 

Since the pion has zero spin, the vertex is extraordi- 
narily simple and takes the form given in Eq. ([2]). Here 
we show how this function is extracted from the ratio 
Eq. dTTJ) by evaluating the correlation functions at large 
Euclidean times. 



First we define the matrix elements as, 



<n| X/1 (o)k(p)> = —= K{p) , 
V 2E *(p) 

Mp)| x t(0)|Q) = J— K® ■ 



Here A w (p) and A w (p) are the couplings of the interpola- 
tor to the pion with momentum p at the sink and source 
respectively. The momentum dependence allows for the 
use of smeared fermion sources and sinks. The bar allows 
for different amounts of smearing at the source and sink. 



By inserting a complete set of energy and momentum states into Eq. (| 13[) . we can show that at large Euclidean 
time, 

& G te>* = ^K®**®*"® ' (18) 

Following the same treatment, one can show that the three-point function Eq. (|15|) at large Euclidean time is 

p-^(p')(*2-tl)p--Ex(p)tl 

lim G a (t2,ti,p',p) = K(p')(n(p')\J a (0)\AP)) ■ (19) 

ti.ta-ti-oo 2 V / E„ (p) En ip>) 



r 



Substituting these expressions into Eq. (fTTj) and using 
Eq. J2]), the ratio R a (p' ,p) simply reduces to 



R a (p',p) 



[p a +P a ']F„(Q 2 ), (20) 



such that the large Euclidean time limits of the ratio R a 
is a direct measure of F 7r (Q 2 ) up to kincmatical factors. 



2. p-meson case 

Following [28| , we define the matrix element of the elec- 
tromagnetic current for p-meson in terms of the covariant 
vertex functions Gi i2 ,3 as in Eqs. (|3|) and ([4]). 



The analogues of the matrix elements in Eq. (fT7|) are 

= A p (p)e M (p,s) 



{n\ X ,(0)\p(p,s)) = 

V 2E p(p) 

(p(p,s)\XuH0)P) = X P (?)et(p,s) .(21) 



The polarisation vectors obey the transversality condi- 
tion 



X/M(P> S ) e t(P: S ) = - 



9nu 



(22) 



because the vector meson current is a conserved current. 



J 



The evaluation of the two- and three-point functions proceeds as for our discussion of the pion. However the 
completeness relation includes a sum over spin-states. Using the transversality condition Eq. (|22p the analogue of 
Eq. (fT8j) becomes 

e -E P (p)t 2 

lim G(t 2 ,p) = > of /-a X p(P) X p(P) e p(P> s ) e t(P, s ) 

to — ►OC * * 



2E p {p) 

e -E p (p)t 2 



2E p (p) 



X p(p) X p(p) 9V 



(23) 



Similarly using Eq. ([3]) we can evaluate the three-point function, 

e --Bp(p')(*2-ti) e -B p (p)t 1 

tuA-^^to^? = L 4E p (p)E p (p>) X " {P ' 



\ p {p')e p (pi,s') e';(p', S ')J Ta °(p' lP )e a (p, S ) \ P {p)e* u {p,s) 
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--E P (p')(*2-ti) e -B P (R)«i 

AE p {p)E p {p>) 



^p(p')^p(P) \9pr 



PuPr 



g<yu 



PaPv 



III: 



(24) 



r 



Inserting the above expressions into the ratio in 
Eq. (|12[) . together with our choice of momentum used 
in the simulations, namely p' = (E p ,p x , 0,0) (E p = 

ri^pj) and p = (m p , 0, 0, 0), it is possible to express 

i?" in terms of the Sachs form factors, 



R u 



G Q (Q 2 ) + ^±^G C (Q 2 ). 



at zero momentum transfer. Since we perform our cal- 
culations at a single, finite value of Q 2 , we will need to 
adjust our results to zero momentum transfer. 

From studies of nucleon properties, it is observed that 
Gm and Gc have similar Q 2 -scaling at small Q 2 [Hj]. In 
the following, we shall assume that this scaling also holds 
for quark contributions to mesons. If Gc(0) = 1, we have 



3m PV /E p m p 



2VE P 



H 13 - H 31 - 



6m py > r E l 



--G Q {Q') 



2v^ 



Px 



2v^ 



=Gm(Q 2 ). 



The individual form factors are isolated as follows: 



Gc(Q 2 ) = (R°n + R° 22 + R° 33 ) , (25) 



3 E p + m p 



Gm(Q 2 ) - ^^(Rls + RlJ, 

Px 



(26) 



Gq (Q 2 ) = "VV^p^p {2R n _ R 2 _ R 3) (27) 

Px 

While we have used the subscript p to denote a vector 
meson, the results are applicable to vector mesons in gen- 
eral, including the K* for example. 



C. Extracting static quantities 

The mean squared charge radius (r 2 ) is obtained from 
the charge form-factor through the following relation, 



^ ^ - 6 w G(Q2) 



Q 2 =0 



(28) 



To calculate the derivative the monopole form is used 

1 \ 



Gc(Q 2 ) 




(29) 



A is referred to as the monopole mass. Inserting this 
form into Eq. (f2"5| and rearranging provides 



(r 2 ) 



1 



() 



Q 2 \G C (Q 2 ) 



(30) 



valid for quantities with Gc{Q 2 = 0) = 1. 

As mentioned in Sec. Ill Al the charge (Eq. ([5])), 
magnetic moment (Eq. ©), and quadrupole moment 
(Eq. (fTU|) ) can be extracted from the Sachs form factors 



Gm(0) 



G m (Q 2 ) 
Gc(Q 2 ) 



(31) 



--G C (Q 2 ), 



Whilst a similar scaling could be used to relate our 
quadrupole form-factor to the quadrupole moment, we 
believe that the form-factor at our small finite Q 2 ( ~ 
0.22GcV) will be of greater phenomcnological interest. 
We note that for a positively charged meson a negative 
value of Gq corresponds to an oblate deformation. 



III. METHOD 

The electromagnetic form factors are obtained using 
the three-point function techniques established by Lein- 
weber, et al. in Refs. [3(J HH, Hil and updated for smeared 
sources in Ref. [3l| . Our quenched gauge fields are gener- 
ated with the 0(a 2 ) mean-field improved Luschcr-Weisz 
plaquette plus rectangle gauge action (35j using the pla- 
quette measure for the mean link. We use an ensemble of 
379 quenched gauge field configurations on 20 3 x 40 lat- 
tices with lattice spacing a = 0.128 fm. The gauge field 
configurations are generated via the Cabibbo-Marinari 
pseudo-heat-bath algorithm [36[ using a parallel algo- 
rithm with appropriate link partitioning [37(. 

We use the fat-link irrelevant clover (FLIC) Dirac op- 
erator [38j which provides a new form of nonperturba- 
tive O(a) improvement |39j |. The improved chiral prop- 
erties of FLIC fcrmions allow efficient access to the light 
quark- mass regime (40| . making them ideal for dynami- 
cal fermion simulations now underway [4l| . For the vec- 
tor current, we an (D(a)-improved FLIC conserved vector 
current [3l|. We use a smeared source at t 2 = 8. Com- 
plete simulation details are described in Ref. (3l| . 

Table [J provides the k- values used in our simulations, 
together with the calculated pseudoscalar and vector me- 
son masses. While we refer to m 2 in our figures and ta- 
bles to infer the quark masses, we note that the critical 
value where the pion mass vanishes is k c1 - = 0.13135. Im- 
portantly the vector mesons remain bound at all quark 
masses considered in this calculation due to finite volume 
effects. That is, the mass of the vector mesons is less than 
the energy of the lowest lying multi-hadron state with the 
appropriate quantum numbers. 

The strange quark mass is chosen to be the third heavi- 
est quark mass. This provides a pseudoscalar mass of 697 
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MeV which compares well with the experimental value of 
(2M 2 K - Ml) 1 ! 2 = 693 MeV motivated by chiral pertur- 
bation theory. Two vector-meson interpolating fields are 
considered, namely q^iq and 97^74(7. Since results for the 
two interpolators agree, we simply present the results 
for the q"fiq interpolator, which displays a significantly 
stronger signal. 

The error analysis of the correlation function ratios 
is performed via a second-order, single-elimination jack- 
knife, with the x 2 P er degree of freedom (xlof) obtained 
via covariance matrix fits. We perform a series of fits 
through the ratios after the current insertion at t = 14. 
By examining the Xdof we are a ^ e to establish a valid 
window through which we may fit in order to extract our 
observables. In all cases, we required a value of Xdof n0 
larger than 1.5. The values of the static quantities quoted 
in this paper on a per quark-sector basis correspond to 
values for single quarks of unit charge. 



TABLE I: Meson masses for the respective values of the hop- 
ping parameter k. 



arriK 



arriK* 



0.12780 
0.12830 
0.12885 
0.12940 
0.12990 
0.13205 
0.13060 
0.13080 



0.5411(10) 
0.5013(11) 
0.4539(11) 
0.4014(12) 
0.3471(15) 
0.3020(19) 
0.2412(42) 
0.1968(52) 



0.4993(11) 
0.4782(11) 
0.4539(11) 
0.4285(11) 
0.4044(12) 
0.3862(13) 
0.3671(19) 
0.3574(16) 



0.7312(30) 
0.7067(36) 
0.6797(46) 
0.6537(49) 
0.6309(56) 
0.6160(64) 
0.6039(71) 
0.5982(80) 



0.7057(27) 
0.6933(40) 
0.6796(46) 
0.6668(47) 
0.6556(50) 
0.6484(52) 
0.6423(54) 
0.6393(56) 



obtained in the previous study reflects finite-volume ef- 
fects attendant with the use of a small spatial volume. 

By comparing the results for the up-quark contribu- 
tions to the 7r and K (p and K*) charge radii, it is pos- 
sible to gain insights into the effect that the presence of 
a heavier strange-quark has on the lighter up-quark in 
pscudoscalar (vector) mesons. Figures [2] and [3] show the 
quark sector contributions to the charge radii ((r 2 )) of 
the pseudoscalar and vector mesons, respectively. The 
quark sector contributions to the charge radii for the 
pscudoscalar and vector meson are recorded in Tables HT1 
and IIII1 From Fig. [5J we find no evidence of environmen- 
tal sensitivity in the light-quark contribution the pseu- 
doscalar mesons. However in the vector sector, Fig. [3J 
we find a consistently broader distribution of up-quark 
charge in the p compared to the up-quark in the K* at 
the smaller quark masses. The broadening of the charge 
distribution in the p is consistent with the hyperfine re- 
pulsion discussed above. The strange quark in the K* 
shows a particularly interesting environment sensitivity. 
While the strange quark mass is held fixed, the distribu- 
tion broadens as the light-quark regime is approached. 
This is consistent with the prediction of enhanced hyper- 
fine repulsion as one of the quarks becomes light. 

The strange neutral pseudoscalar and vector meson 
mean squared charge radii obtained from the weighted 
sum of the quark sector radii are displayed in Fig. 3J For 
the neutral strange mesons, we see a negative value for 
(r 2 ), indicating that the negatively charged d-quark is 
lying further from the centre of mass on average than 
the s. We should expect just such a behaviour for two 
reasons, both stemming from the fact that the s quark is 
considerably heavier than the d: the centre of mass must 
lie closer to the s, and the d-quark will also have a larger 



IV. RESULTS 



Charge radii 



We begin the discussion of our results with the charge 
radii of the vector and pscudoscalar mesons. From the 
quark model we would expect a hyperfine interaction be- 
tween the quark and anti-quark of the form ° q . The 

TTlqTTlq 

interaction is repulsive where the spins are aligned, as 
in the vector mesons, and attractive where the spins are 
anti-aligned, as in the pseudoscalar mesons. In Fig. [1] 
we show the charge radii of the vector and pseudoscalar 
mesons. For comparison the charge radius of the proton 
is also shown. Indeed we find that the charge radii of 
the vector mesons are consistently larger than the pseu- 
doscalar mesons, and in fact similar to the charge radii 
of the proton, even at heavier quark masses. This is con- 
trary to earlier lattice simulations with relatively small 
spatial extent (42j, that have suggested that the ir + , p + 
and proton should have a very similar RMS charge radius 
at larger quark masses. It is possible that the agreement 




0.8 



FIG. 1: Strange and non-strange meson mean squared charge 
radii for charged pseudoscalar and vector mesons. We also 
include for comparison results for the proton taken from 
Ref. [3l[. The ir and p- meson results are centred on the rel- 
evant value of m^., other symbols are offset horizontally for 
clarity. 
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Mean-square charge radius ((r 

2 „.2 



TABLE II 

unit charge in units of fm 
input quark mass. 



2 )) for quarks of 
is given as a measure of the 



m 2 (GeV 2 ) 


II — 




S T.' 


0.6956(26) 


0.216(5) 


0.215(7) 


0.242(7) 


0.5970(26) 


0.225(6) 


0.224(7) 


0.241(7) 


0.4895(24) 


0.240(8) 


0.240(8) 


0.240(8) 


0.3828(23) 


0.256(10) 


0.257(9) 


0.239(9) 


0.2862(25) 


0.274(14) 


0.275(11) 


0.239(10) 


0.2166(27) 


0.287(22) 


0.289(12) 


0.241(11) 


0.1382(48) 


0.304(44) 


0.303(14) 


0.243(13) 


0.0920(49) 


0.287(63) 


0.306(15) 


0.241(13) 
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TABLE III: Mean-square charge radius ((r 2 )) for quarks of 
unit charge in units of fm 2 . 



ml (GeV 2 ) 


Up 


U K ' 


SK* 


0.6956(26) 


0.268(9) 


0.271(11) 


0.309(12) 


0.5970(26) 


0.287(11) 


0.290(13) 


0.311(14) 


0.4895(24) 


0.315(16) 


0.315(16) 


0.315(16) 


0.3828(23) 


0.350(23) 


0.342(20) 


0.321(19) 


0.2862(25) 


0.397(36) 


0.372(26) 


0.331(23) 


0.2166(27) 


0.436(46) 


0.395(30) 


0.339(25) 


0.1382(48) 


0.492(72) 


0.417(35) 


0.353(28) 


0.0920(49) 


0.546(97) 


0.436(41) 


0.360(29) 



Compton wavelength. Of course with exact isospin sym- 
metry in our simulations, the non-strange charge neutral 
mesons have a zero electric charge radius. 

To measure the environmental sensitivity of the light- 
quark sector more precisely, in Figs. [5] and [6] we show 
a fit to the ratio of the light-quark contributions to the 
pseudoscalar and vector-mesons charge radii respectively. 
The difference is striking: for the pseudoscalar case we 
see no environment-dependence at all, whereas in the vec- 
tor case we see that the presence of a strange quark acts 
to heavily suppress the light charge distribution. This 
is the effect one predicts from a quark model, where the 
large mass of the s would act to suppress the hypcrfinc 
repulsion between the quark and anti-quark. It is also 
qualitatively consistent with effective field theory where 
the couplings of the light mesons are suppressed by the 
presence of the strange quark. 

2. Magnetic moments 

In Fig. [7] we present our results for the magnetic mo- 
ments of the vector mesons. At the 5 , C/(3)fl aV our limit, 
where we take the light quark flavours to have the same 
mass as the strange quark, quark model arguments sug- 
gest the magnetic moment for a p + should be -3 times 
the strange magnetic moment of the A (assuming no en- 



FIG. 2: The quark sector contributions to the mean squared 
charge radius of the pseudoscalar mesons. The symbols are 
offset horizontally for clarity. 
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FIG. 3: As for Figure [2] but for vector-mesons. 

vironmental dependence). According to the particle data 
group (43|, the magnetic moment of the A is — 0.613 \x n. 
Therefore we would naively expect a value of 1.84// n for 
the magnetic moment of the p + , which is consistent with 
our findings. 

In an earlier study, Anderson et al. [H| argued that the 
magnetic moment of the p-mcson in natural magnetons 
(otherwise called the g-factor) should be approximately 
2 at large quark masses. Converting our result to natural 
magnetons, we observe in Fig. [S] that our calculation of 
the p-meson g-factor (g p ) is fairly consistent with this. 
At light quark masses, however, we do see some evidence 
of chiral curvature, which would indicate that the linear 
chiral extrapolations of that paper should be considered 
with caution. 

In Fig. [9] we present the quark sector contributions to 
the vector meson magnetic moments, the data is recorded 
in Table. IIV1 Here we observe a similar scenario to that 
observed earlier in the charge radius discussion, namely 
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m 2 (GeV 2 ) 



0.8 



FIG. 4: The mean squared charge radii for the neutral K° 
and K°*. 
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FIG. 5: The ratio of the light quark contributions to the 7r 
and K mean squared charge radius. 



that the u-quark contribution to the K* is consistently 
larger than the contribution from the heavier s-quark. 
We also find that the contribution of the it-quark to the 
magnetic moment of a vector meson is suppressed when 
it is an environment of a heavier s-quark compared to 
when it is in the presence of another light quark. This 
is further supported when we consider the ratio of the 
contributions of a u-quark to the magnetic moments of 
the p and K* mesons, displayed in Fig. [10] This ratio is 
clearly greater than 1 below the S , t/(3)fl av0 ur limit and is 
increasing for decreasing it-quark mass. 

The magnetic moment of the vector meson, like the 
RMS charge radius, shows considerable environment de- 
pendence in the quark sector contributions. The larger 
contribution of a it-quark in a p relative to a K* is consis- 
tent with what we have already observed with the RMS 
charge radius, as follows: since (r 2 ) is larger for the u- 
quark in a p meson than for the it-quark in a K*, the 
effective mass is reciprocally smaller for the u-quark in 
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FIG. 6: As in Fig. [5] but for the vector-mesons. 
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FIG. 7: Charged vector meson magnetic moments. 



a p. This smaller effective mass gives rise in turn to 
a larger magnetic moment. Figure [TU] shows this pat- 
tern. Figure [11] presents our results for the magnetic 
moment of the neutral K*° meson. As the (i-quark be- 
comes lighter than the s we see the magnetic moment 
exhibiting a very linear negative slope. The magnitude 
of the magnetic moment is quite small, but clearly diffcr- 
cntiable from zero everywhere except at the SU (3)fl avour 
limit where symmetry forces it to be exactly zero. 



3. Quadrupole form-factors 

The quadrupole form-factors of the p + and K* + 
mesons are shown in Fig[T2l Wc find that the quadrupole 
form factor is less than zero, indicating that the spatial 
distribution of charge within the p and K* mesons is 
oblate. This is in accord with the findings of Alcxandrou 
et al. [2(| who observed a negative quadrupole moment 
for spin ±1 p- meson states in a density-density analy- 
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FIG. 8: The (/-factor of the p meson. 
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TABLE IV: Magnetic moment for 


quarks of unit char 


je inside 




a vector meson in units of nuclear 


magnetons /xjv- 






(LieV ) 


Up 


U K * 


s K * 




0.6956(26) 


1.71(2) 


1.73(2) 


1.82(3) 




0.5970(26) 


1.77(2) 


1.78(3) 


1.83(3) 




0.4895(24) 


1.84(3) 


1.84(3) 


1.84(3) 




0.3828(23) 


1.94(4) 


1.92(4) 


1.86(3) 




0.2862(25) 


2.04(6) 


1.99(5) 


1.88(4) 




0.2166(27) 


2.11(8) 


2.04(5) 


1.90(4) 




0.1382(48) 


2.20(11) 


2.10(6) 


1.92(5) 




0.0920(49) 


2.25(15) 


2.14(7) 


1.93(5) 
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FIG. 11: Neutral i\*-meson magnetic moment. 
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sis. We note that in a simple quark model, a negative 
FIG. 9: Quark-sector contributions to the vector meson mag- quadrupole form factor requires that the quarks possess 



netic moments. 



an admixture of s- and c?-wavc functions. 




FIG. 10: 

magnetic 
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The ratio of the light-quark contributions to the 
moment of the p and K* . 



TABLE V: The quadrupole form-factor (in fm 2 ) for quarks of 
unit charge inside a vector meson. 



ml (GeV) 2 


u p 


U K ' 


SK* 


0.6956(26) 


-0.0047(4) 


-0.0047(4) 


-0.0048(5) 


0.5970(26) 


-0.0047(4) 


-0.0047(5) 


-0.0048(5) 


0.4895(24) 


-0.0048(5) 


-0.0048(5) 


-0.0048(5) 


0.3828(23) 


-0.0049(7) 


-0.0049(6) 


-0.0048(6) 


0.2862(25) 


-0.0050(9) 


-0.0050(7) 


-0.0048(7) 


0.2166(27) 


-0.0049(12) 


-0.0052(9) 


-0.0047(7) 


0.1382(48) 


-0.0051(19) 


-0.0055(11) 


-0.0046(9) 


0.0920(49) 


-0.0050(27) 


-0.0056(13) 


-0.0046(10) 



The quark sector contributions to the quadrupole 
form- factor are shown in Fig. 1131 The corresponding data 
is contained in Table FVl The flavour independence of the 
results is remarkable. 

We also find that the ratio of the light-quark contri- 
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FIG. 12: Vector meson quadrupole form factors for p + and 
K*+. 
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FIG. 15: Quadrupole form-factor for neutral K* meson. 
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FIG. 13: Quark-sector contributions to the quadrupole form 
factors. 
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FIG. 14: Environment- dependence for light-quark contribu- 
tion to vector meson quadrupole form-factor. 



butions to the quadrupole form factor, shown in Fig. I14| 
is consistent with one within our statistics. In Fig. 1151 
we show the quadrupole form factor of the charge neu- 
tral K*° meson. We find that the quadrupole moment 
of the K*° is non-trivial but just outside the one stan- 
dard deviation level. The chiral trend towards positive 
values reflecting the negative charge of the larger d-quark 
contributions. 

The lattice data for the quark sector contributions to 
the charge form factor is contained in Tables IVII and 
IVIII for the pseudoscalar and vector mesons respectively. 
The magnetic and quadrupole form factors of the vector 
mesons is contained in Tables [Villi and ITXl respectively. 



V. CONCLUSIONS 

We have established a formalism for determining the 
charge, magnetic and quadrupole Sachs form factors of 
vector mesons in lattice QCD. For the first time the elec- 
tric, magnetic, and quadrupole form factors of the light 
vector mesons have been calculated. The electric form 
factor of the pseudoscalar mesons have also been calcu- 
lated. 

With a large lattice volume and high statistics we have 
resolved a clear difference between the charge radii of the 
pseudoscalar and vector mesons. We argue that this is 
consistent with quark model predictions. Furthermore, 
we find significant environmental sensitivity of the light- 
quark contributions to the charge radii of the vector 
mesons. 

We also presented a calculation of the magnetic mo- 
ments of the vector mesons. We found that the magnetic 
moment of the p + was consistent with the quark model 
predication of 1.84 p^ at the 5 , C/(3)fl a vour limit. We de- 
termine that there is also an environmental sensitivity 
in the magnitude of the light-quark contributions to the 
charged vector meson magnetic moments. We argue that 
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this is consistent with the environmental sensitivity in 
the light-quark contributions to the charge vector meson 
charge radii. 

Finally, we have determined that the quadrupole form 
factor for a charged vector meson is negative in quenched 
Lattice QCD. This is consistent with previous calcula- 
tions using density-density analysis. We find that the 
ratio of quadrupole moment to mean square charge ra- 
dius is 1:30, so the deformation is small but statistically 
significant. 

VI. APPENDIX 



TABLE VII: As in Fig. EH but for the vector mesons. 



2 e t t\ 2 

(GeV) 


Up 


UK* 


SK* 


0.6956(26) 


0.795(5) 


0.794(7) 


0.771(7) 


0.5970(26) 


0.784(7) 


0.783(8) 


0.771(8) 


0.4895(24) 


0.769(9) 


0.769(9) 


0.769(9) 


0.3828(23) 


0.750(12) 


0.754(11) 


0.766(10) 


0.2862(25) 


0.727(18) 


0.738(13) 


0.760(12) 


0.2166(27) 


0.708(22) 


0.727(15) 


0.756(13) 


0.1382(48) 


0.683(32) 


0.716(17) 


0.749(15) 


0.0920(49) 


0.660(40) 


0.707(19) 


0.745(15) 



TABLE VI: The quark sector contributions to the charge form 
factor of the pseudoscalar mesons. 



ml (GeV) 2 




UK 


SK 


0.6956(26) 


0.833(4) 


0.835(4) 


0.818(4) 


0.5970(26) 


0.828(4) 


0.830(5) 


0.820(5) 


0.4895(24) 


0.822(5) 


0.822(5) 


0.822(5) 


0.3828(23) 


0.815(6) 


0.813(6) 


0.823(6) 


0.2862(25) 


0.810(8) 


0.804(6) 


0.825(6) 


0.2166(27) 


0.809(12) 


0.798(7) 


0.826(7) 


0.1382(48) 


0.812(22) 


0.792(8) 


0.826(8) 


0.0920(49) 


0.833(30) 


0.791(8) 


0.828(8) 
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